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CLAIMS 

1. A method of quantitative determination of the phase of a radiation 

wave field including the steps of 

(a) producing a representative measure of the rate of change of intensity 
of said radiation wave field over a selected surface extending 
generally across the wave field; 

(b) producing a representative measure of intensity of said radiation 
wave field over said selected surface; 

(c) transforming said measure of rate of change of intensity to produce a 
first integral transform representation and applying to said first 
integral transform representation a first filter corresponding to the 
inversion of a first differential operator reflected in said measure of 
rate of change of intensity to produce a first modified integral 
transform representation; 

(d) applying an inverse of said first integral transform to said first 
modified integral transform representation to produce an 
untransformed representation; 

(e) applying a correction based on said measure of intensity over said 
selected surface to said untransformed representation; 

(f) transforming the corrected untransformed representation to produce a 
second integral transform representation and applying to said second 
integral transform representation a second filter corresponding to the 
inversion of a second differential operator reflected in the corrected 
untransformed representation to produce a second modified integral 
transform representation; 

(g) applying an inverse of said second integral transform to said second 
modified integral transform representation to produce a measure of 
phase of said radiation wave field across said selected plane. 
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2. A method as claimed in claim 1 wherein said first and second integral 
transforms are produced using a Fourier transform. 

3. A method as claimed in claim 2 wherein said Fourier transform is a 
Fast Fourier transform. 

5 4. A method as claimed in any one of claims 1 to 3 wherein said first 

and second differential operators are second order differential 
operators. 

5. A method as claimed in any one of claims 1 to 4 wherein said first 

filter is substantially the same as said second filter. 
io 6. A method as claimed in any one of claims 1 to 5 wherein said first 

filter includes selectively suppressing first higher frequencies of the 
first integral transform representation. 

7. A method as claimed in any one of claims 1 to 5 wherein at least one 
of said first and second filters includes a correction for noise in said 

15 representative measure of intensity. 

8. A method as claimed in any one of claims 1 to 7 including the step of 
producing said representative measures of intensity and rate of 
change of intensity over said selected surface by producing 
representative measurements corresponding to intensity over at least 

20 two spaced apart surfaces extending across the wave field. 

9. A method as claimed in claim 8 wherein said selected surface is 
between two of said spaced apart surfaces. 

10. A method as claimed in claim 8 wherein said selected surface is one 
of said spaced apart surfaces. 

25 11. A method as claimed in any one of claims 8 to 10 including the step 
of directly detecting representative measures of intensity over said 
spaced apart surfaces. 
12. A method as claimed in any one of claims 8 to 10 including the step 
of producing said representative measure of intensity over at least 

30 one of said spaced apart surfaces by imaging that surface. 
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13. 
14. 
5 15. 

16. 

io 

17. 

I5 

18. 

20 19. 



A method as claimed in any one of claims 8 to 12 wherein said 
spaced apart surfaces are substantially parallel. 
A method as claimed in claim 13 wherein said spaced apart surfaces 
are substantially planar. 

A method as claimed in any one of claims 8 to 14 wherein said 
representative measure of rate of change of intensity is produced by 
subtraction of representative measurements of intensity respectively 
made at locations over said spaced apart surfaces. 
A method as claimed in any one of claims 1 to 15 wherein said 
representative measures of intensity and rate of change of intensity 
are obtained by sampling measurements at selected locations over 
said surface. 

A method as claimed in claim 16 wherein, said sampling 
measurements are made at locations defining a regular array over 
said surface. 

A method as claimed in claim 2 or claim 3 wherein said radiation 
wave field propagates in a z-direction of a Cartesian co-ordinate 
system and further including the step of producing an x component 
and a y component of phase separately. 

A method as claimed in claim 18 wherein said first and said second 
filters have a component Q x for producing the x component of phase 

and a component Cl y for producing the y component of phase of the 



form 



(k 2 x +k 2 y )k x 



(k]+k 2 y ) 2 +ak] 



{k] +kpk y 



(kl+k]) 1 +ak] 



where k x k y are the Fourier variables conjugate to x and y; 
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a is a constant determined by noise in the intensity measurements. 

20. A method as claimed in claim 19 including the step of multiplying said 
representative measure of rate of change of intensity by the negative 

5 of average wave number of the radiation before said integral 

transformation. 

21 . A method as claimed in any one of claims 1 to 6 including the step of 
obtaining said representative measure of rate of change of intensity 
by obtaining a first representative measurement over a measurement 

io surface across the wave field for radiation of a first energy and 

obtaining a second representative measurement over said 
measurement surface for radiation of a second different energy. 

22. A method as claimed in claim 2 or claim 3 wherein at least one of said 
first filter and said second filter include a correction for aberrations in 

15 said representative measures of intensity and rate of change of 

intensity by including at least one component dependent on the 
aberration, coefficients of a system producing the representative 
measures. 

23. A method as claimed in claim 22 wherein said radiation wave field 
20 propagates in a z-direction of a Cartesian co-ordinate system and 

further including the step of producing an x component and a y 
component of phase separately. 

24. A method as claimed in claim 23 wherein said first and said second 
filters have a component Q x for producing the x component of phase 

25 and a component Q y for producing the y component of phase both of 

the form 



30 



wher k Xl ky are the Fourier variables conjugate to x and y; 
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X is the average wave length of the radiation; 

/aberrated(x,y) is the aberrated intensity measured at defocus distance 
&, 

A mn are the aberration coefficients which characterize the imperfect 
5 imaging system. 

25. A computer program to execute the steps of any one of claims 1 to 
24. 

26. A computer program stored on computer readable storage media 
including means to execute the steps of any one. of claims 1 to 24. 

10 27. An apparatus for quantitative determination of the phase of a 

radiation wave field including 
(a) means to produce a representative measure of the rate of change of 

intensity of said radiation wave field over a selected surface 
extending generally across the wave field; 
is (b) means to produce a representative measure of intensity of said 

radiation wave field over said selected surface; 
(c) processing means to sequentially 

(I) transform said measure of rate of change of intensity to 
produce a first integral transform representation; 
20 (II) apply to said first integral transform representation a first filter 

corresponding to the inversion of a first differential operator 
reflected in said measure of rate of change of intensity to 
produce a first modified integral transform representation; 

(III) apply an inverse of said first integral transform to said first 
25 modified integral transform representation to produce an 

untransformed representation; 

(IV) apply a correction based on said measure of intensity over said 
selected surface to said untransformed representation; 

(V) transform the corrected untransformed representation to 
30 produce a second integral transform representation; 
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(VI) apply to said second integral transform representation a 
second filter corresponding to the inversion of a second 
differential operator reflected in the corrected untransformed 
representation to produce a second modified integral transform 

5 representation; and 

(VII) apply an inverse of said second integral transform to said 
second modified integral transform representation to produce a 
measure of phase of said radiation wave field across said 
selected plane. 

io 28. An apparatus as claimed in claim 27 wherein said first and second 

integral transforms are produced using a Fourier transform. 

29. An apparatus as claimed in claim 28 wherein said Fourier transform is 
a Fast Fourier transform. 

30. An apparatus as claimed in any one of claims 27 to 29 wherein said 
is first and second differential operators are second order differential 

operators. 

31. An apparatus as claimed in any one of claims 27 to 30 wherein said 
first filter is substantially the same as said second filter. 

32. An apparatus as claimed in any one of claims 27 to 31 wherein said 
20 first filter includes selectively suppressing first higher frequencies of 

the first integral transform representation. 

33. An apparatus as claimed in any one of claims 27 to 31 wherein at 
least one of said first and second filters includes a correction for noise 
in said representative measure of intensity. 

25 34. An apparatus as claimed in any one of claims 27 to 33 including 

means to produce representative measurements corresponding to 
intensity over at least two spaced apart surfaces extending across the 
wave field. 

35. An apparatus as claimed in claim 34 wherein said selected surface is 

30 between two of said spaced apart surfaces. 
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36. An apparatus as claimed in claim 34 wherein said selected surface is 
one of said spaced apart surfaces. 

37. An apparatus as claimed in any one of claims 34 to 36 including 
detector means positioned to directly detect representative measures 

5 of intensity over said spaced apart surfaces. 

38. An apparatus as claimed in any one of claims 34 to 36 including 
detector means to produce said representative measure of intensity 
over at least one of said spaced apart surfaces and imaging means to 
image that surface onto the detector. 

10 39. An apparatus as claimed in any one of claims 34 to 38 wherein said 

spaced apart surfaces are substantially parallel. 

40. An apparatus as claimed in any one of claims 34 to 38 wherein said 
spaced apart surfaces are substantially planar. 

41. An apparatus as claimed in any one of claims 34 to 40 wherein said 
is means to produce said representative measure of rate of change of 

intensity subtracts representative measurements of intensity 
respectively made at locations over said spaced apart surfaces. 

42. An apparatus as claimed in any one of claims 27 to 41 wherein said 
means to produce a representative measure of intensity and said 

20 means to produce a representative measure of rate of change of 

intensity sample at selected locations over said surface. 

43. An apparatus as claimed in claim 42 wherein said samples are made 
at locations defining a regular array over said surface. 

44. An apparatus as claimed in claim 28 or claim 29 wherein said 
25 radiation wave field propagates in a z-direction of a Cartesian co- 
ordinate system and processing means produces an x component 
and a y component of phase separately. 

45. An apparatus as claimed in claim 44 wherein said processing means 
applies said first and said second filters have a component Q x for 

30 producing the x component of phase and a component Ci y for 

producing the y component of phase of the form 
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(k:+k 2 y )k x 



{k]+k 2 y ) 2 +ak] 



(k*+k 2 y )k y 



(k 2 x +k 2 y ) 2 +ak] 



46. 

10 

47. 

15 

48. 

20 

25 49. 



where k X( k y are the Fourier variables conjugate to x and y; 

a is a constant determent by noise in the intensity measurements. 

An apparatus as claimed in claim 37 wherein said representative 

measure of rate of change of intensity is multiplied by the negative of 

the average wave number of the radiation before said integral 

transformation. 

An apparatus as claimed in any one of claims 27 to 31 wherein said 
representative measure of rate of change of intensity is produced by 
obtaining a first representative measurement over a measurement 
surface across the wave field for radiation of a first energy and 
obtaining a second representative measurement over said 
measurement surface for radiation of a second different energy. 
An apparatus as claimed in claim 28 or claim 29 wherein at least one 
of said first filter and said second filter include a correction for 
aberrations in said representative measures of intensity and rate of 
change of intensity by including at least one component dependent 
on the aberration, coefficients of a system producing the 
representative measures. 

An apparatus as claimed in claim 48 wherein said radiation wave field 
propagates in a z-direction of a Cartesian co-ordinate system and 
wherein an x component and a y component of phase are produced 
separately. 
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50. An apparatus as claimed in claim 49 wherein said first and said 

second filters have a component Q x for producing the x component of 

phase and a component Q v for producing the y component of phase 

both of the form 



1 



where k x ,k y are the Fourier variables conjugate to x and y; 
X is the average wave length of the radiation; 
jo /aberrated(x,y) is the aberrated intensity measured at defocus distance 

&, 

>A mn are the aberration coefficients which characterize the imperfect 
imaging system. 

51 . A method of imaging an object including the steps of 

is (a) exposing the object to a radiation wave field from a source; 

(b) producing a representative measure of the rate of change of 
intensity over a selected surface extending generally across the 
wave field on the side of the object remote from incident 
radiation; 

20 (c) producing a representative measure of intensity of said radiation 

wave field over said selected surface; 

(d) transforming said measure of rate of change of intensity to 
produce a first integral transform representation and applying to 
said first integral transform representation a first filter 

25 corresponding to the inversion of a first differential operator 

reflected in said measure of rate of change of intensity to 
produce a first modified integral transform repr sentation; 

(e) applying an inverse of said first integral transform to said first 
modified integral transform repres ntation to produce an 

30 untransformed representation; 
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(f) applying a correction based on said measure of intensity over 
said selected surface to said untransformed representation; 

(g) transforming the corrected untransformed representation to 
produce a second integral transform representation and 
applying to said second integral transform representation a 
second filter corresponding to the inversion of a second 
differential operator reflected in the corrected untransformed 
representation to produce a second modified integral transform 
representation; 

(h) applying an inverse of said second integral transform to said 
second modified integral transform representation to produce a 
measure of phase of said radiation wave field across said 
selected plane. 

52. A method as claimed in claim 51 including the step of producing said 
15 representative measures of intensity and rate of change of intensity 

over said selected surface by producing representative 
measurements corresponding to intensity over at least two spaced 
apart surfaces extending across the wave field. 

53. A method as claimed in claim 52 wherein said selected surface is 
20 between two of said spaced apart surfaces. 

54. A method as claimed in claim 52 wherein said selected surface is one 
of said spaced apart surfaces. 

55. A method as claimed in any one of claims 52 to 54 wherein said 
spaced apart surfaces are substantially parallel. 

25 56. A method as claimed in any one of claims 52 to 54 wherein said 

spaced apart surfaces are substantially planar. . 
57. A method as claimed in any one of claims 52 to 56 wherein said 

representative measure of rate of change of intensity is produced by 
subtraction of representative measurements of intensity respectively 

30 made at locations over said spaced apart surfaces. 
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58. A method as claimed in claim 51 including the step of producing said 
measures of intensity and rate of change of intensity over said 
selected surface by producing first representative measurements 
corresponding to intensity over a first surface extending across the 

5 wave field, changing the distance between said source and said 

object, and producing second representative measurements 
corresponding to intensity over said first surface for the changed 
distance between said object and said source. 

59. A method as claimed in claim 58 wherein said selected surfaces is 
io one of said spaced apart surfaces. 

60. A method as claimed in any one of claims 51 to 57 including the step 
of directly detecting representative measures of intensity over said 
spaced apart surfaces. 

61. A method as claimed in any one of claims 51 to 60 wherein said 
15 selected surface is spaced apart from said object in the direction of 

propagation of said radiation. 

62. A method as claimed in any one of claims 51. to 61 wherein said 
source is substantially a point source. 

63. A method as claimed in any one of claims 51Jo 62 wherein said first 
20 and second integral transforms are produced using a Fourier 

transform. 

64. A method as claimed in any one of claims 51 to 63 wherein said 
Fourier transform is a Fast Fourier transform. 

65. An apparatus for imaging an object including : 

25 (a) a source to irradiate the object with a radiation wave field; 

(b) means to produce a representative measure of the rate of 
change of intensity of said radiation wave field over a selected 
surface extending generally across the wave field. 

(c) means to produce a representative measure of intensity of said 
30 radiation wave field over said selected surface; 

(d) processing means to sequentially 
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(i) transform said measure of rate of change of intensity to 
produce a first integral transform representation; 

(ii) apply to said first integral transform representation a first 
filter corresponding to the inversion of a first differential 
operator reflected in said measure of rate of change of 
intensity to produce a first modified integral transform 
representation; 

(iii) apply an inverse of said first integral transform to said first 
modified integral transform representation to produce an 
untransformed representation; 

(iv) apply a correction based on said measure of intensity over 
said selected surface to said untransformed 
representation; 

(v) transform the corrected untransformed representation to 
produce a second integral transform representation; 

(vi) apply to said second integral transform representation a 
second filter corresponding to the inversion of a second 
differential operator reflected in the corrected 
untransformed representation to_ produce a second 
modified integral transform representation; and 

(vii) apply an inverse of said second integral transform to said 
second modified integral transform representation to 
produce a measure of phase of said radiation wave field 
across said selected plane. 

An apparatus as claimed in claim 65 including means to produce 
representative measurements corresponding to intensity over at least 
two spaced apart surfaces extending across the wave field. 
An apparatus as claimed in claim 66 wherein said selected surface is 
between two of said spaced apart surfaces. 

An apparatus as claimed in claim 66 wherein said selected surface is 
one of said spaced apart surfaces. 
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69. An apparatus as claimed in any one of claims 66 to 68 including 
detector means positioned to directly detect representative measures 
of intensity over said spaced apart surfaces. 

70. An apparatus as claimed in any one of claims 66 to 68 including 
detector means to produce said representative measure of intensity 
over at least one of said spaced apart surfaces and imaging means to 
image that surface onto the detector. 

71. An apparatus as claimed in any one of claims 66 to 70 wherein said 
spaced apart surfaces are substantially parallel. 

72. An apparatus as claimed in any one of claims 66 to 71 wherein said 
representative measure of rate of change of intensity is produced by 
subtraction of representative measurements of intensity respectively 
made at locations over said spaced apart surfaces. 

73. An apparatus as claimed in claim 65 including means to produce said 
measures of intensity and rate of change of intensity overs said 
selected surface by producing first representative measurements 
corresponding to intensity over a first surface extending across the 
wave field; means to change the distance between said source and 
said object, and means to produce ^second representative 
measurements corresponding to intensity over said first surface for 
the changed distance between said object and said source. 

74. An apparatus as claimed in claim 73 wherein said selected surface is 
one of said spaced apart surfaces. 

75. An apparatus as claimed in any one of claims 65 to 72 including 
means to directly detecting representative measures of intensity over 
said spaced apart surfaces. 

76. An apparatus as claimed in any one of claims 65 to 75 wherein said 
selected surface is spaced apart from said object in the direction of 
propagation of said radiation. 

77. An apparatus as claimed in any one of claims 65 to 76 wherein said 
source is substantially a point source. 
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78. An apparatus as claimed in any one of claims 65 to 77 wherein said 
first and second integral transforms are produced using a Fourier 
transform. 

79. An apparatus as claimed in claim 78 wherein said Fourier transform is 
5 a Fast Fourier transform. 

80. A method of phase amplitude imaging including the steps of 

(a) irradiating an object with a radiation wave field; 

(b) focussing radiation from the object through an imaging system 
to an imaging surface extending across the wave field 

1D propagating from the object; 

(c) producing a first representative measure of intensity 
distribution of radiation over said imaging surface at a first 
focus of the imaging system; 

(d) introducing a change in focus of the image on said imaging 
15 surface through the imaging system; 

(e) producing a second representative of measure intensity 
distribution over said imaging surface; and 

(f) using said first and second representative measures to 
produce a representative measure^ of intensity and a 

20 representative measure of rate of change of intensity in the 

direction of radiation propagation over a selected surface 
extending across the wave field; 

(g) transforming said measure of rate of change of intensity to 
produce a first integral transform representation and applying 
to said first integral transform representation a first filter 
corresponding to the inversion of a first differential operator 
reflected in said measure of rate of change of intensity to 
produce a first modified integral transform representation; 

(h) applying an inverse of said first integral transform to said first 
30 modified integral transform representation to produce an 

untransformed representation; 



25 
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(i) applying a correction based on said measure of intensity over 

said selected surface to said untransformed representation; 
(j) transforming the corrected untransformed representation to 
produce to said second integral transform representation a 
second integral transform representation and applying a 
second filter corresponding to the inversion of a second 
differential operator reflected in the corrected untransformed 
representation to produce a second modified integral transform 
representation; 

(k) applying an inverse of said second integral transform to said 
second modified integral transform representation to produce a 
measure of phase of said radiation wave field, across said 
selected plane. 

81. A method as claimed in claim 80 wherein said radiation wave field 
has a numerical aperture smaller than the numerical aperture of said 
imaging system. 

82. A method as claimed in claim 80 or claim 81 said first focus of the 
imaging system produces an infocus image at the imaging surface 
and said second focus of the imaging system produces a slightly 

20 defocused image at the imaging surface. 

83. A method as claimed in any one of claims 80 to 82 wherein said 
imaging surface is substantially planar. 

84. A method as claimed in any one of claims 80 to 83 wherein the 
imaging surface is an intensity detector. 

25 85. A method as claimed in any one of claims 80 to 84 wherein said 

imaging surface is said selected surface. 

86. A method as claimed in any one of claims 80 to 85 wherein said 
integral transform is a Fourier transform. 

87. A method as claimed in claim 86 wherein said Fourier transform is a 
30 Fast Fourier transform. 



15 
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A method as claimed in any one of claims 80 to 87 wherein said 
representative measure of rate of change of intensity is produced by 
subtraction of said first and second representative measurements of 
intensity. 

89. A method as claimed in any one of claims 80 to 88 wherein said 

representative measures of intensity and rate of change of intensity 
are obtained by sampling measurements at selected location over 
said imaging surface. 

A method as claimed in claim 89 wherein said sampling 
measurements are made at locations defining a regular array over 
said imaging surface. 

91. A method as claimed in claim 87 wherein said radiation wave field 
propagates in a z-direction of a cartesian co-ordinate system and 
further including the step of producing an x component and a y 

15 component of phase separately. 

92. A method as claimed in claim 91 wherein said first and said second 
filters have a component Q x for producing the x component of phase 
and a component Cl y for producing the y component of phase of the 



90 

10 



25 



form 

-2 . Jr 2> 



(k 2 x+ k 2 y )k x 

20 i2x = 



(k 2 x +k 2 y )k y 
(k 2 x +k 2 y ) 2 +ak. 



where k x .k y are the Fourier variables conjugate to x and y; 
a is a constant determined by noise in the intensity measurements. 
93. A method as claimed in claim 92 including the step of multiplying said 

representative measure of rate of change of intensity by the negative 
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of the average wave number of the radiation before said integral 
transformation. 

94. An apparatus for phase amplitude imaging of an object including 

a radiation wave field source to irradiate said object; 
5 an imaging system to focus radiation from said object to an 

imaging surface extending across the wave field propagating from the 
object; 

means to produce a representative measure of radiation 
intensity over said imaging surface; 
io said imaging system including selectively operable means to 

adjust said focus of said radiation to said imaging surface to at least 
at a first focus and a second focus; 

processing means to: 
(i) produce a representative measure of intensity and a 
is representative measure of rate of change of intensity in the 

direction of radiation propagation over a selected surface 
extending across the wave field from representative 
measures of radiation intensity over said image surface at 
said first focus and said second focus ; 
20 (ii) transform said measure of rate of change of intensity to 

produce a first integral transform representation; 

(iii) apply to said first integral transform representation a first filter 
corresponding to the inversion of a first differential operator 
reflected in said measure of rate of change of intensity to 

25 produce a first modified integral transform representation; 

(iv) apply an inverse of said first integral transform to said first 
modified integral transform representation to produce an 
untransformed representation; 

(v) apply a correction based on said measure of intensity over 
30 said selected surface to said untransformed representation; 
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(vi) transform the corrected untransformed representation to 
produce a second integral transform representation; 

(vii) apply to said second integral transform representation a 
second filter corresponding to the inversion of a second 

5 differential operator reflected in the corrected untransformed 

representation to produce a second modified integral 
transform representation; and 

(viii) apply an inverse of said second integral transform to said 
second modified integral transform representation to produce 

io a measure of phase of said radiation wave field across said 

selected plane. 

95. An apparatus as claimed in claim 94 wherein said radiation wave field 

has a numerical aperture smaller than the numerical aperture of said 
imaging system. 

15 96. An apparatus as claimed in claim 94 or claim 95 said first focus of the 

imaging system produces an infocus image at the imaging surface 
and said second focus of the imaging system produces a slightly 
defocused image at the imaging surface. 

97. An apparatus as claimed in any one of claims 94 to 96 wherein said 
20 imaging surface is substantially planar. 

98. An apparatus as claimed in any one of claims 94 to 97 wherein the 
imaging surface is an intensity detector. 

99. An apparatus as claimed in any one of claims 94 to 98 wherein said 
imaging surface is said selected surface. 

25 100. An apparatus as claimed in any one of claims 94 to 99 wherein said 
integral transform is a Fourier transform. 

101. An apparatus as claimed in claim 100 wherein said Fourier transform 
is a Fast Fourier transform. 

102. An apparatus as claimed in any one of claims 94 to 101 wherein said 
30 representative measure of rate of change of intensity is produced by 
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subtraction of said first and second representative measurements of 
intensity. 

103. An apparatus as claimed in any one of claims 94 to 102 wherein said 
representative measures of intensity and rate of change of intensity 

5 are obtained by sampling measurements at selected location over 

said imaging surface. 

104. ' An apparatus as claimed in claim 103 wherein said sampling 

measurements are made at locations defining a regular array over 
said imaging surface. 

io 105. An apparatus as claimed in claim 101 wherein said radiation wave 
field propagates in a z-direction of a cartesian co-ordinate system and 
further including the step of producing an x component and a y 
component of phase separately. 
106. An apparatus as claimed in claim 105 wherein said first and said 
second filters have a component Q x for producing the x component 
of phase and a component n y for producing the y component of 
phase of the form 



= {k]+ki)k x 



15 



20 



(k 2 x +k 2 ) 2 +ak 2 



(k 2 x +kpk y 
y (k 2 x +k 2 y ?+akl 



where k Xl k y are the Fourier variables conjugate to x and y; 
a is a constant determent by noise in the intensity measurements. 
25 107. A method as claimed in claim 106 including the step of multiplying 
said representative measure of rate of change of intensity by the 
negative of the average wave number of the radiation before said 
integral transformation. 
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108. A computer program for processing a representative measure of the 
rate of change of intensity of a radiation wave field over a selected 
surface extending generally across the wave field and a 
representative measure of intensity of said radiation wave field over 
5 said selected surface, said program including: 

(a) code for transforming said measure of rate of change of intensity to 
produce a first integral transform representation and applying to said 
first integral transform representation a first filter corresponding to the 
inversion of a first differential operator reflected in said measure of 

io rate of change of intensity to produce a first modified integral 

transform representation; 

(b) code for applying an inverse of said first integral transform to said first 
modified integral transform representation to produce an 
untransformed representation; 

15 ( c ) code for applying a correction based on said measure of intensity 

over said selected surface to said untransformed representation; 

(d) code for transforming the corrected untransformed representation to 
produce a second integral transform representation and applying to 
said second integral transform representation a second filter 

20 corresponding to the inversion of a second differential operator 

reflected in the corrected untransformed representation to produce a 
second modified integral transform representation; 

(e) code for applying an inverse of said second integral transform to said 
second modified integral transform representation. 

25 109. A computer program stored on computer readable storage media for 
processing a representative measure of the rate of change of 
intensity of a radiation wave field over a selected surface extending 
generally across the wave field and a representative measure of 
intensity of said radiation wave field over said selected surface, said 

30 program including: 
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(a) code for transforming said measure of rate of change of intensity to 
produce a first integral transform representation and applying to said 
first integral transform representation a first filter corresponding to the 
inversion of a first differential operator reflected in said measure of 
rate of change of intensity to produce a first modified integral 
transform representation; 

(b) code for applying an inverse of said first integral transform to said first 
modified integral transform representation to produce an 
untransformed representation; 

(c) code for applying a correction based on said measure of intensity 
over said selected surface to said untransformed representation; 

(d) code for transforming the corrected untransformed representation to 
produce a second integral transform representation and applying to 
said second integral transform representation a second filter 
corresponding to the inversion of a second differential operator 
reflected in the corrected untransformed representation to produce a 
second modified integral transform representation; 

(e) code for applying an inverse of said second integral transform to said 
second modified integral transform representation. 



